Edwards-like statistical mechanical description of the parking lot model for vibrated 

granular materials. 
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We apply the statistical mechanical approach based on the "flat" measure proposed by Edwards 
and coworkers to the parking lot model, a model that reproduces the main features of the phe- 
nomenology of vibrated granular materials. We first build the flat measure for the case of vanishingly 
small tapping strength and then generalize the approach to finite tapping strengths by introducing 
a new "thermodynamic" parameter, the available volume for particle insertion, in addition to the 
particle density. This description is able to take into account the various memory effects observed in 
vibrated granular media. Although not exact, the approach gives a good description of the behavior 
of the parking-lot model in the regime of slow compaction. 

PACS numbers: 05.70.Ln,45.70.Cc 
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I. INTRODUCTION 

Granular media are a-thermal, out-of-equilibrium sys- 
tems that it would be useful to describe within a statis- 
tical mechanical framework. A given macro-state of such 
a system characterized by a fixed density of grains (con- 
sider for simplicity a packing of monodisperse spherical 
particles) is very likely to be associated with an expo- 
nentially large number of micro-states or particle con- 
figurations. How the packing was prepared (by pour- 
ing, shaking, shearing, etc..) may influence its proper- 
ties and change the way the associated particle config- 
urations are sampled when repeating over the same ex- 
perimental protocol. However, the simplest hypothesis, 
put forward by Edwards and his coworkers (ll y, 0, 
is that all micro-states characterized by a given av- 
erage density are equiprobable. With this "flat mea- 
sure" , on can build a statistical mechanical framework 
in which entropy, i.e., the logarithm of the number 
of micro-states, is the relevant thermodynamic poten- 
tial. This approach has recently been the focus of an 
intense research activity, in connection with a series 
of experimentsperformed on weakly vibrated granular 
materials!^ 11,0 and with a theoretical description 
of out-of-equilibrium glassy systems based on the concept 
of effective temperature [13 El IH El El ■ 

In the past few years, the Edwards' hypothesis has 
been tested on many models, virtually all of them b eing 
lattice models with some kind of "tapping" kinetics ^a, 

S III El IM El El IM 111 IM 1^ 

ISU 133 . l33 | . In the absence of experimental tests of this 
approach (see, however,^), such theoretical studies are 
expected to better circumscribe the conditions of validity 
of the statistical mechanical description. (Presumably, 
only "approximate validity" can be expected since, aside 
from specific mean-field modelsjTolll4|. such a simplified 
description of out-of equilibrium situations in terms of a 
small number of "thermodynamic" parameters is unlikely 



to be exact.) 

In this article, we consider an Edwards-like sta- 
tistical mechanical description approach for the one- 
dimensional model of random adsorption-desorption of 
hard particles .3_4, 35, 36] , also known as the parking- lot- 
modeljl|. This latter is a microscopic, off-lattice model 
that mimics many features of the compaction of a vi- 
brated column of grains. Besides a qualitative descrip- 
tion of the phenomen olog y of weakly tapped granular 
mediaH H IS El El ElllJ, the interest of the model 
is that exact analytical results can be derived or, when 
not possible, very accurate numerical data can be ob- 
tained from computer simulations. In the next section, 
we briefly introduce the parking-lot model and we dis- 
cuss its connection to vibrated granular materials. In 
section HI, we consider the limit of vanishingly small (but 
non-zero) tapping intensity; we construct for this case an 
Edwards-like description based on a flat measure over all 
"blocked" states and we compare the resulting predic- 
tions to the exact behavior. In the following section, we 
generalize the study to the case of a finite tapping inten- 
sity: we consider what appears to be the simplest, yet 
compatible with known experimental observations, gen- 
eralization of the Edwards' formalism. Finally, we discuss 
the merits and limitations of the approach. 



II. THE MODEL AND ITS CONNECTION TO 
VIBRATED GRANULAR MATERIALS 

The parking-lot model is a one-dimensional random 
adsorption-desorption process of hard rods on a line. 
Hards rods of length cr are deposited at random posi- 
tions on a line at rate fc+ and are effectively inserted 
if they do not overlap with pre-deposited rods; other- 
wise they are rejected. In addition, all deposited par- 
ticles can desorb, i.e., be ejected from the line at ran- 
dom with a rate Time is measured in units of l/fc+, 
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length in units of a, and the model depends on one con- 
trol parameter K — k^/k_. When no desorption is 
present — 0), the model reduces to the purely ir- 
reversible one-dimensional random sequential adsorption 
(RSA) process 4^ 42c i4J|i also known as the car park- 
ing problem, and all the properties of the system as a 
function of time are known exactly |3 EM 119 • In addi- 
tion, for 1 / K non strictly equal to zero, the competition 
of mechanisms between adsorption and desorption allows 
the system to reach a steady state that is nothing but an 
equilibrium fluid of hard rods at constant activity 1/ K: 
there too, all properties are known exactly. 

The densification kinetics of the parking-lot model at 
constant K is described by 



dp 
dt 



K 



^ ' K 



(1) 



where p{t) the density of hard rods on the line at time 
t (recall that cr = 1) and is the fraction of the line 
that is available at time t for inserting a new particle, i.e., 
the probability associated with finding an interval free of 
particles (a "gap") of length at least 1. The quantities 
p and $ can be calculated from the 1— gap distribution 
function G{h, t) which is the density of gaps of length h 
at time t via a number of "sum rules" : 



Pit) = 
1 - p{t) = 

m = 



dhG{h,t), 
dhhGih,t), 



dh{h - l)G{h,t). 



(2) 
(3) 
(4) 



The evolution with time of the 1— gap distribution func- 
tion can itself be described by a kinetic equation that 
involves 2— gap distribution functions, and so on 37]. Ex- 
cept for the two above mentioned limits (RSA when 
k- = 0, equilibrium when t — s- -l-oo), the infinite 
hierarchy of coupled equations cannot be solved an- 
alytically and one must resort to approximate treat- 
ments and computer simulations, as described in previous 
articles 37, 38, 39, 40]. 

First introduced in the context of protein adsorption at 
liquid-solid interfaces [3 HE ES, the random adsorption- 
desorption model has recently been applied to the de- 
scription of weakly vibrated granular materials Q 0, 113, 
HE HE EE ^] ■ The connection between the parking- lot 
model and these latter is made by regarding the particles 
on the line as an average layer of grains in the vibrated 
column. Time measures the number of taps whose effect 
is to eject particles from the layer; ejection is followed 
by the arrival at random of particles in the layer, which 
mimics the gravity-driven relaxation step in the experi- 
ment. Considering that the main influence of the inten- 
sity of the tapping is to determine the average number 
of particles ejected at each tap, (this number being an 
increasing function of intensity), leads to associate 1/K 
with the tapping strength. A two-dimensional version 



of the model with some polydispersity of the particles 
would clearly be more realistic, but one does not expect 
this to change the qualitative features of the model |43|. 
A more serious caveat is the absence of explicit account 
of the mechanical stability of the particle packings: sta- 
bility is only implicitly described by the fact that the 
particles are blocked on the line between two successive 
desorption events. 

Despite its drastic simplification of the situation en- 
countered in vibrated granular materials, the parking lot 
model reproduces at a qualitative level most of the rel- 
evant phenomenology: (i) for large rate K correspond- 
ing to weak tapping intensity, compaction proceeds very 
slowly and can be effectively described by an inverse 
logarithm of timejSJ, |33, |33, yg; (ii) stronger tapping 
leads to faster initial compaction but to less effective 
asymptotic packingj^^; (hi) the slow densification ki- 
netics leads to irreversibility effects and to the obser- 
vation of two curves for the packing density as a func- 
tion of tapping intensitv|37j one essentially reversible 
and one irreversible depending on the experimental pro- 
tocol chosen,22j; (iv) the power spectrum of the den- 
sity fluctuations near the steady state is distinctly non- 
Lorenztian and displays a power-law regime at interme- 
diate frequencies 6, 7, 39l l4fl I4H : (v) non-trivial memory 
effects are observed when changing abruptly the tapping 
intensitv p9lLi4Q] . 

Note that, as recently studied for a one-dimensional 
model with tapping dynamics |3Cl|. the parking lot model 
could also be used, via the introduction of two kinds of 
particles, to describe the segregation phenomena with the 
so-called Brazil Nut jiTl lisj and Reverse Brazil Nut^E 
^1 efl:ects. 



III. LIMIT OF VANISHINGLY SMALL 
TAPPING INTENSITY: K ^ 00 

In this limit, ejection of one particle from the line is 
followed by an infinite number of insertion trials until 
one, or seldom, two new particles are added. The stable 
or "blocked" configurations are thus those for which no 
more particle insertions are possible (recall that once suc- 
cessfully inserted particles cannot move on the line), i.e., 
all configurations of non-overlapping rods such that the 
available line fraction $ is zero, or, equivalently, such that 
all gaps between neighboring particles are smaller than 
a particle size (here taken as unity). Edwards' prescrip- 
tion for constructing a statistical-mechanical description 
of this system is then to consider that all such "blocked" 
configurations at a fixed density p are equiprobable. 

Consider a line of length L with N particles. With 
periodic boundary conditions, this system has also 
N gaps between neighboring particles. Denoting by 
/ii, /i2, . . . , /lAf the lengths of these gaps, the total num- 
ber of "blocked" configurations is given by the config- 
urational integral calculated under the constraint that 
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hi < 1 for i ~ 1, .., iV, namely, 



Z{L, ^) ^'^O S{L-N-J2 h.), (5) 



which by using the integral representation of the delta- 
function, can be rewritten as 



dhie 



'Zhi 



where C denotes a closed contour. Integrating over the 
h[s yields 



Z{L,N) = j dzexp ( L 



z(l - p) + pin 



1 — exp(— z) 



(7) 



where p = N/ L. In the macroscopic limit where N — > 
oo, L — > oo, with fixed p, the above expression can be 
evaluated though a saddle-point method, which gives 

Z(L,iV)~exp(is(p)) (8) 

where s(p), the entropy density, is expressed as 



s{p) = (l-p)2; + pln 



(9) 



where z = z{p) is the solution of the saddle-point equa- 
tion 



1 



1 



z 1 — e 



(10) 



In Edwards' language, ^ = (^M^ ) com- 

pactivity (up to a trivial constant) j5l|. In an equilibrium 
system of hard rods, i.e., a flat measure without the con- 
straint that all gaps have a length smaller than 1, z(p) 
would simply be equal to P/^kgT) ~ p/(l — p) where P 
is the pressure. 

By Legendre transforming the entropy S{L,N), one 
obtains a new potential 



Y{N, z) = -S(L, N) + zL = N 



-In 



1-e- 



(11) 



such that 



djY) 
dz 



N 



{L) and from which one can obtain 



the fluctuations of the system size. 



(L^) - {Lf = 



d'^iY) 



dz^ 



= -N 



N 



dz 



(12) 
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By combining the above expression with Eq. H1U|) . one 
derives the fluctuations of the density. 



L{{p')-{pr)^i^{{L')~{Lr) 



= p 



z2 (l-e-^)2 



(13) 



In the Edwards-like ensemble, one can also calculate 
the gap distribution functions. The 1-gap distribution 
function G{h; p), which gives the density of gaps of length 
h, is obtained from 



GEd{h; p) 



1 



LZ{L,N) 



(6) E 



Et-ti n dhAs\L-N-h- J2 
=1-^0 Jo J \ ,^1,,^, J 



(14) 



By using the same method as before one finds that 



-, — zh 



GEdih;p) = ri-e- 

b 



for h <1, 
for /i > 1. 



(15) 



where z is the solution of Eq. (|l()|l . It is easy to check that 
the above expression satisfies the two sum rules, Eqs. ^ 
and 

The higher-order gap distribution functions arc ob- 
tained along the same lines, and they satisfy a factoriza- 
tion property analogous to that found for an equilibrium 
system of hard rods, i.e., 



GEd{h, h'; p) = GEd{h; p)GEd{h'; p) 



(16) 



GEd{h,h\h"-p) = GEd{h-p)GEd{h'-p)GEd{h"-p), 

(17) 

etc. 

The 1-gap distribution function G{h\ p) is directly re- 
lated to the nearest-neighbor pair distribution function 
that represents the probability of finding two neighboring 
particles whose centers are separated by a distance l-\-h. 
It is also possible to calculate the full pair distribution 
function g{r] p) via a method which closely follows that 
developed for the equilibrium system of hard rods'^. 
The steps of the calculation are detailed in Appendix A, 
and the final result reads 



9Ed[ 



(-^ p) E ^^7^^ ( E ^^(-i)'^^(^ - - - fc) 



m— 1 



p ^-^ (m — 1) 

^ m=l ^ ' 



1 - e- 



-z{r — rn) 



(18) 



where 9{x) is the Heaviside step function and z is the 
solution of Eq. H1U|I . For comparison, we give the equi- 
librium pair distribution function|52||. 

9e,{r;p)^-}^ ^?(r-m)^-— ^ 



e(-(T^)('^-™)) (19) 

1 -P/ 
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FIG. 1: Density fluctuations of the parking lot model when 
K — > +00 for densities above the RSA jamming limit. Simu- 
lation data are shown as the full line, the Edwards approxi- 
mation corresponds to the dashed curve, and the equilibrium 
result to the dotted curve. The inset displays the density fluc- 
tuations at equilibrium and in the Edwards approximation for 
a larger range of densities. 



We can now compare the above results derived under 
the condition of equiprobabihty of the "blocked" config- 
urations with the exact ones obtained either analytically 
or numerically. When K ^ oo, analytical results are 
available in the two limits, t — 0^ , which corresponds to 
the purely irreversible RSA process at the jamming limit 
where no more particles can be inserted |44j. and t — > -|-oo 
which corresponds to a close- packed state with p = 1. 

For the RSA at the jamming limit, closed- form expres- 
sions have been derived for the saturation density |4Jj|, 

PJL= dt exp ^-2 '^'^—^ — ) - 0.74759 . . . , 

(20) 

for the density fluctuations 45j, for the gap distribution 
functions !46l| and for the pair distribution function[45l| . 
(The expressions are given in Appendix B.) When com- 
paring to the Edwards-like results at the same density, 
Pjl, one finds qualitative differences. Most notably, (i) 
the exact 1-gap distribution function displays a logarith- 
mic divergence at contact between particles {h — + 0^), 

G{h;pjL)^-e-^''Hh) (21) 

where 7 is the Euler constant, (ii) the exact multi-gap 
distribution functions do not reduce to products of the 
1-gap functions, (iii) the exact pair distribution function 



FIG. 2: Entropy density versus particule density for a hard- 
rod system: Edwards entropy for the parking lot model when 
A" ^ 00 (full curve) and equilibrium entropy (dashed curve). 

has a super-exponential decay at large distances, 

5(r,pji)-l^-i^, r^oo (22) 

r(r) 

all features that are missed by the fiat-measure expres- 
sions, since 

GEd{h = 0;pjL)=PjLj^^^j^, (23) 

1 - z{pjL) 

gEd{r,PjL)-l--e'^''''^^\ r^oo (24) 

and the multi-gap functions satisfy a factorization prop- 
erty, Eqs. 

Quantitatively, one can also see differences, e.g., in 
the density fluctuations, the exact result at jamming be- 
ing L((p2) - (p)2) ~ 0.038 to be compared to i((p^) - 
{p)^)Ed ^ 0.028. 

Such observations, that generalize to an off-lattice 
model the results obtained by de Smedt et al. |3ll . |35| for 
random and cooperative sequential adsorption models on 
a one-dimensional lattice, are in fact to be expected: it 
has been shown that the RSA process generates config- 
urations of hard particles that are sampled from a prob- 
ability distribution that is "biased" when compared to a 
fiat measure (53ii54] . 

As the system further evolves with time at vanishingly 
small tapping intensity, compaction takes place beyond 
the RSA saturation density, and the difference between 
the actual properties of the system and those predicted 
by the Edwards-fiat measure diminishes. 

Figurencompares the variation with p (for p > pjl) of 
the density fluctuations of the parking lot model (simu- 
lation results) with the flat-measure result, Eq. H13|) . and 
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FIG. 3: Parking lot model when K ^ +00: log-linear plot of 
the 1-gap distribution function versus h for several densities: 
from top to bottom p = 0.75, 0.85, 0.89, 0.91. Comparison of 
simulation data (wavy line) and Edwards approximation (full 
curve) . 



FIG. 4: Pair distribution function in the parking lot model 
when K —* +00 a,t p = 0.77: simulation data (wavy curve), 
Edwards approximation (full curve) and equilibrium (dashed 
curve) . 



the equilibrium curve, L ((p'^) — {p)'^)^^ — p{^~ pY- The 
flat-measure prediction is in fair agreement with the sim- 
ulation data, especially for intermediate densities, but it 
reaches too rapidly the equilibrium curve that (slightly) 
overestimates the exact result. (Note that simulating the 
system becomes very time consuming as compaction goes 
on and, in practice, we cannot go beyond a density of 
about 0.93; the infinite-time limit should of course be 
p = 1). However the Edwards description improves upon 
the equilibrium curve and gives the proper shape of the 
density dependence with an inflection around p ~ 0.87. 
The closeness of the flat measure result and of the equilib- 
rium one (the former being of course with the constraint 
that no gaps have a length larger than 1) at high density 
is illustrated in Fig. [21 where we show the entropy density 
versus particle density. 

In Figure Owe have displayed on a logarithmic- linear 
plot the predicted and computed 1-gap distribution func- 
tions for four different densities. The Edwards approx- 
imation is too small at small h for p = 0.75, which is 
reminiscent of the missing the logarithm divergence at 
the RSA jamming limit pjL — 0.747 (see above). Oth- 
erwise the agreement with the simulation data is very 
good (recall that this is a log- linear plot). For /i > 1, the 
Edwards approximation is exact since GEd{h > 1) = 0. 

Finally, we have plotted in Figures ^ the pair dis- 
tribution functions (simulation data, Edwards approxi- 
mation, equilibrium curve) for three different densities, 
p = 0.77, p = 0.82 and p = 0.92. The (constrained) flat 
measure is an improvement upon equilibrium curve, but 
even at p = 0.82 (Fig. [S]), it somewhat overestimates the 
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FIG. 5: Same as Fig.Hfor p = 0.82 



oscillations at large distances. For p — 0.92 (Fig. O, the 
difference between the three curves is barely visible. Note 
that for 1 < r < 2, the pair distribution function is equal 
to the nearest-neighbor distribution function, hence to 
the 1-gap distribution function shown in Fig. |31 
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FIG. 6: Same as Fig.Hfor p = 0.92 



IV. FINITE TAPPING INTENSITY: K FINITE 

When the tappmg intensity, i.e., is finite, one 

must modify the definition of the stable or "blocked" 
states. Indeed, for finite the configuration of hard 
rods obtained after a desorption-adsorption event are no 
longer characterized by $ — 0, and the gaps between 
particles can be larger than one. One also knows that 
the tapping intensity 1/K is not the proper "thermody- 
namic" parameter to add to the density in order to char- 
acterize the macro-state of the system in an Edwards- like 
statistical mechanical approach: the memory effect ob- 
served experimentally [3 and reproduced by the present 
parking lot model (see above) implies that the system can 
be found in states characterized by the same density p 
and the same tapping intensity 1/K that however evolve 
differently under further tapping with the same intensity, 
l/K] as illustrated in Fig.[7| the density may increase in 
one case and decrease in another. 

If the tapping intensity is not an appropriate thermo- 
dynamic parameter, a natural choice for a two-parameter 
statistical mechanical description appears to be $, the 
available line fraction, that one can use in conjunction 
with the density p. A non-zero (f> generalizes to a finite 
tapping intensity, the prescription uses for vanishingly 
small intensity namely $ = 0, and $ is also directly rel- 
evant for describing the compaction kinetics, as shown 
by Eq. We thus consider a statistical mechanical 

ensemble in which all configurations of non-overlapping 
hard rods characterized by fixed values of p and $ are 
equally probable. 

Denoting by A the total length available for insertion 
for a particle center {A = the configurational inte- 
gral with the constraints of fixed A, fixed system size L, 



FIG. 7: Memory effect in the parking lot model. The full 
curve corresponds to a process at constant K — 2000 whereas 
the dashed curve shows the kinetics when K is switched from 
2000 to 500 at ts = 1000. The points A and B corresponds 
to states with equal density, equal value of K — 500, but 
different further evolution. The inset shows the phenomenon 
with a larger scale (upper curve: constant K = 500, lower 
curve: constant K — 2000). 



and fixed number of particles N is given by 



Z{L,N,A) 



-'[■■■[ n dh,6 f L - - V /i, ) 

•^0 -^0 ^=l V . = 1 / 

s(^A-J29{h,-l)ih,-l)^ , (25) 



which can be rewritten as before as 



Z{L, N,A)^ / dz dy 
Jc Jc 



exp L 



z{l - p) + + pin 



z + y{l — cxp(-z)) 
z(z + y) 



(26) 



where C and C denote two closed contours. In the 
macroscopic limit, N —>■ oo, L oo, A oo with p 
and $ fixed, one can again use a saddle point method to 
evaluate the integrals, which leads to 

Z{L,N,A)^eMLs{p,^)) (27) 

where s{p, $) is expressed as 

'z + y{l -exp(-z))\ 



s{p, $) = (1 - p)z + y^ + pin 



z{z + y) 



(28) 



7 



with z = z{p^ $) and y = y{p, $) solutions of the two 
coupled equations 



1-p 



1 



1 + ye~ 



P 



1 

z z + y z + y{l — e~^)' 
1 



z + y z + y{l — e ^) 



(29) 
(30) 



z = {^) N A '^^^^ again be considered as the inverse of the 
"compactivity" , but an additional intensive parameter, 
y= is needed. 

A double Legendre transformed potential Y{N,z,y) 
can be introduced as 



Y{N, z, y) = -Ls{p, $) + + yA 

'z + y{l -exp(-z)) 



= N 



- In 



z{z + y) 



(31) 



Then, ^ 



(L) and ^ 



tuations in L and A are given by 



N,z 



{A), and the flue- 



{Lf 
{A') (Af 
{LA)-{L){A) 



dz^ 
d^{Y) 



d{l/p) 



N,v 



dy^ 
d^{Y) 



-N 



dz 
9($/p) 



N.z 



N 



= -N 



dzdy 
d{l/p) 



dy 
d{^/p) 



dz 



dy 



(32) 



(33) 



(34) 



By using the saddle-point equations, Eqs. (|29|l and 
(|5n)l. one arrives at the following expression for the fluc- 
tuations of /9, 



3 / 1 I 1 l + i2 + z + y)ye- 

^ U2 (z + y)2 ^z + y{l-e-)y 



(35) 



Similar expressions are obtained for the fluctuations 
of $ and the cross fluctuations of p and but are not 
shown here. 

The gap distribution functions can also derived by fol- 
lowing the same method as in the previous section. This 
leads to 



z{z + y) 



, ^ , _ , „ for /i< 1, 

GE,ih;p)^l ^ + y(l-e 

p^ilpd_^-i^-h+yih^l)) fo^ ^^^^ 



z + y {1 — e ^) 



(36) 



whereas the multi-gap distribution functions satisfy 
the factorization property, e.g., G Ed{h,h' ] p,<^) = 
GEd{h; p,^)GEd{h'; p,^)- Notice that the 1-gap distri- 
bution function is a piecewise continuous function that 
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FIG. 8: Density fluctuations as a function of p for K — 500. 
The dotted curve corresponds to the equilibrium density fluc- 
tuations, the dot-dashed curve to the RSA result, the dashed 
curve to the 2-parameter Edwards prediction, and the wavy 
line to the simulation data. The inset displays the same curves 
for a larger range of densities. 
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FIG. 9: Same as Fig.Hfor K = 5000. 



obeys the exact sum rules, Eqs. (0- (@J. It is also worth 
pointing out that the results of section III {K 00) can 
be recovered by taking the limit y 00 in the above 
equations. Finally, the pair distribution function can be 
derived along the same lines as shown before and in Ap- 
pendix A, but the calculation is too tedious and not suf- 
ficiently insightful to be presented here. 
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V. CONCLUSION 




FIG. 10: Parking lot model when K = 500: log-linear 
plot of the 1-gap distribution function versus h for several 
densities: from top to bottom in the middle of the figure 
p = 0.682,0,748,0.790,0.83. Comparison of simulation data 
(wavy line) and Edwards approximation (full curve). For 
p = 0.83, there is virtually no difference between the two 
curves. 



A comparison between the 2-parameter Edwards flat 
measure and the simulation data is shown in Figs.|Hland|3 
for the density fluctuations with K — 500 and K — 5000, 
respectively. We have also plotted the equilibrium curve, 
^{{p'^)eq ~ i£^q) = ^ p)^ i ^^^^1 thc l-dimcnslonal 
RSA curve pSj up to pjL- The 2-parameter flat mea- 
sure predictions are good but not perfect. They display 
the proper non-trivial shape of the p-dependence, in con- 
tradistinction to the equilibrium curve, but there is a 
significant underestimation of the fluctuations in the den- 
sity range around the RSA jamming limit. Note that at 
high density the flat measure is in very good agreement 
with the simulation data (at least for K = 500) and does 
not merge too rapidly with the equilibrium curve as seen 
above for the case K ^ oo. 

The 1-gap distribution functions shown on a log-linear 
plot in Fig. ^1 illustrate also the overall good agree- 
ment between the 2-parameter flat-measure predictions 
and the simulation data. The Edwards approximation 
captures the change of the slope of the 1-gap distribu- 
tion function that occurs for /i > 1; but the curvature 
seen in the simulation data for h slightly larger than 1 is 
not correctly reproduced by the 2-parameter flat-measure 
which predicts an exponential decay with a factor equal 
to 2 -I- y (see Eq. (I36|l ). Again, this discrepancy is larger 
for densities around of the RSA jamming limit (the two 
intermediate sets of curves in Fig. I10|l . 



In this work we have applied the statistical mechan- 
ical approach based on the "flat" measure proposed by 
Edwards and coworkers to the out-of-equilibrium situa- 
tion obtained in the parking-lot model. This latter is 
a microscopic off-lattice model that reproduces the ma- 
jor features of the phenomenology of vibrated granular 
materials. In the statistical mechanical description, a 
macro-state of the system is characterized by fixed val- 
ues of two macroscopic quantities, the particle density p 
and the available line fraction $ (or rather fixed values 
of 3 extensive parameters, the number of particles, the 
system size and the total length available for insertion 
of particles), and all configurations of non-overlapping 
particles with fixed p and are taken as equiprobable. 

We show that such an approach misses some of the 
qualitative signatures of the limiting case of a purely irre- 
versible adsorption process (RSA) at the jamming limit. 
However, at higher densities, i.e., in the slow (logarith- 
mic) compaction regime, it gives a good, yet not per- 
fect, quantitative description of many observable quan- 
tities (fluctuation, distribution functions). The choice of 
$ as an additional "thermodynamic" parameter is able 
to account for situations, encountered in various memory 
effects, in which macrostates characterized by the same 
density and the same tapping intensity can nonetheless 
be different. 

The fact that a "thermodynamic" approach gives a 
good description of a model of vibrated granular media 
is promising. In the one-dimensional model studied here, 
the generalized "equation of state" associated with the 
flat measure can be analytically derived so that one can 
make theoretical predictions concerning, e.g., the density 
fluctuations or the structure of the conflgurations. How- 
ever, one still faces the task of predicting the state of 
the system for a given preparation, i.e., for a given time 
(number of taps) and a given protocol for the tapping 
intensity. We are presently working on this problem. 



APPENDIX A: PAIR DISTRIBUTION 
FUNCTION IN THE EDWARDS ENSEMBLE 
{K oo) AND AT EQUILIBRIUM 

We introduce the probability density '^m{0 of flnd- 
ing two given particles at a relative distance ^, such that 
there is exactly m — 1 particles between them; this func- 
tion can be expressed in terms of the partition function 
as 

ZiC,m-l)ZiL-^,N-m) 
Z{L,N) 

where Z{L, N) can be calculated either with the Edwards 
measure or at equilibrium. 

For large TV, one can use the asymptotic expression 
of the partition function Z{L,N) = e^(^-^)0(z)^. At 
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equilibrium (j){z) = 1/z whereas with the Edwards mea- 
sure, one gets (j){z) = (1 — e^~^^)/z and z given by 
Eq. 1)10(1 . Therefore, the probability density is equal to 

vI/„(C) = Z(e,m-l)^e(-«-'"» (A2) 

and, formally, the pair distribution function can be ex- 
pressed as 



by|4a|4 



_ J2/q°° dttkitfe^-"'^ for h<l, 



for /i > 1, 



(Bl) 



with 



^ oo 



(A3) 



For hard rods, the partition function Z{(_, to — 1) is differ- 
ent from for ^ > to. At equilibrium, one has for ^ > m 



k{t) = exp 



1 - e- 



-du 



(B2) 



. /o " 

The pair distribution function can be expressed in a 
closed form by using the Laplace transform g{s,t) = 

p+co 

/ dle^''''g{l + l,t). At the jamming limit, one has'455 

^0 



(m - 1)! 



(A4) 



and z = p/(l — p), which gives for the pair distribution 
function l53| 



1 °° 

9eq{r,p) =- ^ 



(r — TO)(r — to) 



m— 1 



■ exp 



(A5) 



1 










POO 


- 2 






Jo 










Jo 





dti 



k{ti)k{ti+s) 
k{s) 



k{h)k{ti+s) fc(i2)fc(i2 + S) 



Ms) 
B(s,t3) 



Ms) 
(B3) 



when r > 1, where 9{x) is the Heaviside function. 

In the Edwards ensemble, the expression for Z{^, m — 
1), hence for 4'm{i') involves to -I- 1 contributions. For 
instance, ^i{r) is given by 



^'i(r) = {e{r - 1) -6'(r - 2)) 



1-e- 



z{r-l) 



(A6) 



More generally 

m / 



Inserting Eq. l|Tf)l in Eq. l|X3)l yields Eq. (jTHJ 



(A7) 



APPENDIX B: RSA EXPRESSIONS AT THE 
JAMMING LIMIT 



with 



Bis,t) = 



s + t 



1 - e-(^+*) 
(.s + t)2 



(B4) 



The expression for the density fluctuations follows from 
the above equation by taking the limit s ^ of the 
expression p (l -(- 2p {g{s, pjL)e~^ — 7)) 0|- This gives 



PJL [ -l + 2pji-- J dhk'ih) 



dt2k^{t2) / dt3e-''k-\h) 



0.038. 



1 1 - e^*-' 



(B5) 



For the one-dimensional RSA process, the 1-gap dis- At the same density, pjl, the equilibrium value is 
tribution function at the jamming limit {t 00) is given L ((p^) — (p)^)^ = 0.0476. 
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